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ABSTRACT 



Aims. By means of chemical evolution models for ellipticals, spirals and irregular galaxies, we aim at investigating the physical 
meaning and the redshift evolution of the mass-metallicity relation as well as how this relation is connected with galaxy morphology. 
Methods. Our models distinguish among different morphological types through the use of different infall, outflow, and star formation 
prescriptions. We assume that galaxy morphologies do not change with cosmic time. We present a method to account for a spread in 
(~ i , the epochs of galaxy formation and to refine the galactic mass grid. To do that, we extract the formation times randomly and assume 

^ |. an age dispersion A,. We compare our predictions to observational results obtained for galaxies between redshifts 0.07 and 3.5. 

I ' Results. We reproduce the mass-metallicity (MZ) relation mainly by means of an increasing efficiency of star formation with mass 

O ' in galaxies of all morphological types, without any need to invokegalactic outflows favoring the loss of metals in the less massive 

galaxies. Our predictions can help constraining the slope and the zero point of the observed local MZ relation, both affected by 
^ . uncertainties related to the use of different metallicity calibrations. We show how, by considering the MZ, the O/H vs star formation 

^ I ' rate (SFR), and the SFR vs galactic mass diagrams at various redshifts, it is possible to constrain the morphology of the galaxies 

producing these relations. Our results indicate that the galaxies observed at z = 3.5 should be mainly proto-ellipticals, whereas at 
z = 2.2 the observed galaxies consist of a morphological mix of proto-spirals and proto-ellipticals. At lower redshifts, the observed MZ 
^ ■ relation is well reproduced by considering both spirals and irregulars. Galaxies with different star formation histories may overlap in 

' the MZ diagram, but measures of abundance ratios such as [0/Fe] can help to break this degeneracy. Predictions for the MZ relations 

QQ ' for other elements (C, N, Mg, Si, Fe) are also presented, with largest dispersions predicted for elements produced in considerable 

amounts by Type la SNe, owing to the long lifetimes of their progenitors. 

Key words, galaxies: abundances; ISM: abundances; galaxies:high-redshift; galaxies: evolution. 

o 

1. Introduction low gravitational potential wells of the former (Larson 1974; 

S< ' , , , Dekel & Silk 1986; Tremonti et al. 2004, De Lucia et al. 2004; 

> ■ ye^ars, several observational projects have been Kobayashi, Springel & White 2007, Finlator & Dave 2008). 

devoted to the study of the relatiori between galactic stellar mass Alternatively, the dilution caused by infall can act in a way sim- 
and metalhcity in star forming galaxies. In particular, the study ^.^^^^^^ ^^^j^^^^^ ^^^^ ^^^^^^^ ^ .^^^^^ 

H of the redshift evolution of the mass-metallicity (hereinafter MZ) ^.^^^ j^^^^. ^^^^ j^^^^^ (Dalcanton et al. 2004). Another 
relation has provided us with crucial information on the cosmic ibiUty is that low mass galaxies are less evolved than large 

evolution of star formation, as well as on the tempora evolution j^^^^^ ^^.^ .^^^^^^ ^^^^ j j^^j^^ ^^^^ ^^^^^^ 

of the chemical properties of the stellar populations Various ob- ^^j,^ ^j^^^^. ^^^^^ ^ ^^^^^ .^^^^^^ ^^^^^ ^.^ ^^^^^.^^^ 
servational studies of the MZ relation have outhned that across g^j^hing quickly their ISM to solar or over-solar metallicities, 
a large redshift range, i.e. between z=3.5 and z=0, in the high- ^^^^ j^^j^^^ characterized by lower star formation efficien- 
mass, high-metalhcity part of the MZ plot, galaxies tend to form ^-^^ (• ^ formation rates per unit mass of gas) have sub-solar 
a plateau whereas at smaller masses the MZ plot shows an in- ^^.^^^.^^^^^ metallicities. This interpretation is supported by sev- 
crease of metalhcity with mass. Moreover, both the zero point ^^^j e^^lution models (Lequeux et al. 1979, Matteucci 

and the slope of the MZ relation are a function of redshift, as ^994), by cosmological N-body simulations (Brooks et al. 2007; 
shown in a recent paper by Maiolino et al. (2008). Therefore, ^ouhcine et al. 2008, Tassis et al. 2008) and by hydrodynam- 
the metalhcity evolution is stronger for low mass galaxies than .^^^ simulations (Tissera et al. 2005, De Rossi et al. 2007,). A 
toi high mass galaxies. . . , . . , thirdinterpretationof the MZ relation is linked to the initial mass 

Various theoretical explanations of the MZ relation have beeii ^^^^^.^^ ^ ^2007) showed how the MZ relation can 

proposed so far One interpretation is linked to starburst-induced explained by a higher upper mass cutoff in the initial mass 
galactic outflows, more efficient in expeUing metal-enriched gas ^^^^^.^^ ^j^P^ ^^^^ ^^^^.^^ j^^.^^ 

in dwarf galaxies than i n giant galaxies, owing to the shal- However, none of these theoretical studies have investigated how 
Send offprint requests to: F Calura *e MZ relation is connected with galaxy morphology, as well as 
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the role played by the star formation histories of galiixies of dif- 
ferent morphological types. 

In this paper, our aim is to investigate these aspects, as well 
as to understand the physical meaning of the MZ relation for 
different galaxies. To perform this task, we use detailed chemi- 
cal evolution models for elUptical, spiral and irregular galaxies. 
These models have proven to well reproduce the main chemi- 
cal properties of elliptical (e.g. Pipino & Matteucci 2004), spi- 
rals (Chiappini et al. 2001, 2003; Cescutti et al. 2007) and ir- 
regular galaxies (e.g. Lanfranchi & Matteucci, 2003), as well as 
other observational constraints, such as the gas-to-Ught ratios, 
the supernova rates as a function of morphological type (Calura 
& Matteucci 2006a) and the cosmic evolution of the luminosity 
density (Calura & Matteucci 2003; Calura, Matteucci & Menci 
2004). 

This paper is organized as follows. In Section 2, we present 
our models and our methods to investigate the MZ relation. 
In section 3, we describe the set of observational data used 
in this work and we present our results. Finally, in Section 4 
we draw some conclusions. Throughout the paper, we adopt a 
concordance A-Cold Dark Matter cosmology, characterized by 
0„ = 0.3, Qa = 0.7 andh = 0.7. 

2. Chemical Evolution Models 

In this paper, we use chemical evolution models for elliptical, 
spiral and irregular galaxies that have already been tested against 
various local observational constraints. In general, in chemical 
evolution models, the time-evolution of the fractional mass of 
the element i in the gas within a galaxy, G„ is described by the 
basic equation: 

Gi = -m^iit) + Ri(t) + (Gdinf - (Gdout (1) 

where G,(?) = Mg(t)Xi(t)/M,ot is the gas mass in the form of 
an element / normalized to the total baryonic mass M,ot and 
G(f) = Mg(t)/Mtot is the total fractional mass of gas present in 
the galaxy at the time t. The quantity X,(f) = Gi{t)/G{t) is the 
abundance by mass of an element i, with the summation over aU 
elements in the gas mixture being equal to unity. if/(t) is the star 
formation rate (SFR), namely the fractional amount of gas turn- 
ing into stars per unit time. Ri(t) is the returned fraction of matter 
in the form of an element i (both newly produced and already 
present in the stars) that the stars restore into the ISM through 
stellar winds. Type la and Type II supernova explosions. (G,),„/ 
and (G,)o«/ describe the possible infall of external gas and the 
possible presence of outflows, respectively. 
The nucleosynthesis prescriptions are common to all models. 
For massive stars and Type la SNe, we adopt the empirical yields 
suggested by Francois et al. (2004), which are based on Woosley 
& Weaver (1995) for massive stars and on the Type la SNe 
yields of Iwamoto et al. (1999). Fow low and intermediate mass 
stars, the adopted prescriptions are the ones by van den Hoeck 
& Groenewegen (1997). The Type la SN rate computation is 
based on the single-degenerate model and follows the Matteucci 
& Recchi (2001) prescriptions. 
The SFR tl/{t) is a Schmidt (1959) law expressed as: 

m = vG\t) (2) 

The quantity v is the efiiciency of star formation, namely the 
inverse of the typical time-scale for star formation, and is ex- 
pressed in Gyr 



Unless otherwise stated, the rate of gas infall, for a given element 
i, is defined as: 



with Cinf being a suitable constant, tuned in order to reproduce 
the present-day stellar mass, and r the infall timescale. 
In all our models, the instantaneous recycling approximation 
(IRA) is relaxed. This means that the chemical abundances are 
computed by taking into account the stellar hfetimes. 
In this work, we assume that interstellar abundances are not af- 
fected by dust depletion. This assumption is motivated by the 
fact that the observational abundances considered here are mea- 
sured in Hii regions, where dust grains are destroyed by intense 
the UV radiation fields generated by massive stars (Okada et al. 
2008). Quantitative estimates of the effects of dust on metalUcity 
measures in Hu regions indicate that the presence of dust grains 
does not introduce large errors in the global metallicity indica- 
tors being the uncertainty in the metallicity due to dust efl^ects 
in any case < 0.2 dex (Shields & Kennicutt 1996). Furthermore, 
local depletion measurements in diff'erent environments indicate 
that O should not be considered as a strongly refractory element 
(Jensen et al. 2005) . This provides further support to our as- 
sumption of neglecting dust depletion effects in MZ studies. 
We assume that galaxy morphology does not change with red- 
shift. When discussing difl^erent galaxy types at high redshift, it 
may seem inappropriate to use the classification in terms of ellip- 
ticals, spirals and irregulars, as we observe them today. For this 
reason, when discussing our results at any redshift, whenever 
we use to the expression "eUipticals", "spirals" or "irregulars", 
we refer to the high-redshift counterparts of each morphological 
type, or more appropriately to proto-eUipticals, proto-spirals and 
proto-irregulars. 

2.1. Elliptical galaxies 

For the chemical evolution of elliptical galaxies, we adopt 
the one- zone model of Pipino & Matteucci (2004) where we 
address the reader for further details. Here we assume that all 
spheroids, i.e. elliptical galaxies, spiral bulges and halos are all 
included in the same category, described in this section. Here 
we recall the main assumptions: ellipticals form by means of a 
rapid infall of pristine gas, which triggers an intense starburst. 
Star formation is assumed to halt when the energy of the ISM, 
heated by stellar winds and supernova (SN) explosions, exceeds 
the binding energy of the gas (Dekel & Silk 1986). At this time 
a galactic wind occurs, sweeping away almost all of the residual 
gas present at that time. After the SF has stopped, the galactic 
wind is maintained by Type la SNe which continue exploding 
until the present time, and its duration depends on the balance 
between this heating source and the gas cooling. The binding 
energy of the gas is strongly influenced by assumptions con- 
ceming the presence and distribution of dark matter (Matteucci 
1992). For the model adopted here, a diffuse {Re/Rd-O.l, where 
Re is the effective radius of the galaxy and Rj is the radius of the 
dark matter core) but massive (Mdark/Mium = 10) dark halo has 
been assumed (see Bertin et al. 1992). 

The outflow rate {Gi)oui is of the same order of magnitude 
(in general within a factor of 2) of the value taken by the SFR 
just before the galactic wind (Pipino et al. 2005). This result is in 
agreement with the indications from the observations of starburst 
galaxies (e.g. Heckman 2002). For further details, see Pipino et 
al. (2005). 
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Table 1. Elliptical models: parameters 



Baryonic Mass 




V 


Tinf 


IMF 




(kpc) 




(Gyr) 




10'" 


1 


3 


0.5 


Salpeter 


10" 


3 


10 


0.4 


Salpeter 


10'2 


10 


25 


0.2 


Salpeter 



We assume that the efficiency of star formation is higher in 
more massive objects which evolve faster than less massive 
ones (inverse-wind scenario, Matteucci, 1994, otherwise called 
"downsizing"). This implies that the galactic wind develops af- 
ter a timescale varying with the galactic mass, with more mas- 
sive galaxies experiencing earlier outflows. This mechanism al- 
lows us to well reproduce the observed increase of [Mg/Fe] 
with galactic mass in ellipticals (see Matteucci 1994; Pipino & 
Matteucci 2004). It is also worth noting that this wind mech- 
anism is the opposite of what proposed originally by Larson 
(1974) to explain the MZ-relation in ellipticals. In fact, in 
Larson's model the galactic wind develops later in more mas- 
sive systems due to an assumed constant efficiency of star forma- 
tion with galactic mass. Matteucci (1994) obtained the "inverse- 
wind" scenario simply by assuming an increasing efficiency of 
star formation with mass and such an assumption can preserve 
the MZ-relation as well. 

A Salpeter (1955) IMF constant in space and time is adopted, 
with lower limit O.IM0 and upper limit lOOM©. 

The choice of such an IMF for ellipticals and SO galaxies 
assures that several observational constraints such as the average 
stellar abundances and the color-magnitude diagram (see Pipino 
& Matteucci, 2004) are well reproduced, as well as the metal 
content in clusters of galaxies (see Renzini 2004; Calura, 
Matteucci & Tozzi 2007). We consider three elliptical galaxy 
models of three different luminous masses: M;„„, = 10'°, 10" 
and IO'^Mq. The infall is assumed to occurr on an extremely 
short timescale (< 0.5 Gyr). In the Table 1 we list the initial 
baryonic mass, the effective radius, the star formation effi- 
ciency V, the infall timescale t,„/ and the IMF adopted for 
each galaxy. In Fig. [1] we show the predicted time evolution 
of the star formation rates. Type la. Type II SNe rates and 
interstellar metallicity, represented by the O abundance (in units 
of 12-i-log(0/H) for the three elliptical galaxy models used in 
this paper. 



2.2. Spiral galaxies 

To model spiral galaxies, we use a single-infall model designed 
to describe the evolution of the thin disc of the Milky Way galaxy 
(Chiappini et al. 2001). By using this model, we make the im- 
plicit assumption that the baryonic mass of any spiral galaxy is 
dominated by a thin disc of stars and gas in analogy with the 
Milky Way. The model used for spiral discs is a multi-zone one. 
The disc is approximated by several independent rings, 2 kpc 
wide, without exchange of matter between them. The timescale 
for the disc formation is assumed to increase with the galac- 
tocentric distance, thus producing an "inside-out" scenario for 
the disc formation (Matteucci & Francois 1989; Chiappini et al. 
2001;CescuttietaL 2007). 

In this work, we use three spiral models representing galax- 
ies of different masses. One model is designed to reproduce the 
main features of the MlOl spiral galaxy (see Chiappini et al. 
2003), representing the most massive spiral disc. The MW disc 




time/Gyr 



time/Gyr 



Fig. 1. From top left corner, clockwise: SFR, interstellar O abun- 
dance (in units log(0/H)-Hl2), Type II SNR, and Type la SNR vs 
time for the three elliptical galaxy models used in this paper The 
solid, dotted and dashed lines are the predictions for the models 
with total baryonic mass IO^'Mq, IO^Mq, and IO'^Mq, respec- 
tively. 



is used to represent a spiral galaxy of intermediate mass. Finally, 
a model for a low-mass spiral has been developed, characterized 
by a star formation efficiency and a mass surface density lower 
than the ones of the Milky Way (see Tab. 2). 
Each model is characterized by a particular infall law and a star 
formation efficiency. Our main assumption here is that larger 
discs evolve faster than smaller ones (Boissier et al. 2001), in 
analogy with ellipticals and in agreement with observations. For 
the mass surface density of each spiral disc, we assume an expo- 
nential profile l,tot(R) = £0 exp(-/?//?rf). The SFR surface den- 
sity is a Kennicutt (1998) law expressed as: 



v<.(f) 



(4) 



with k - 1 .5. As for ellipticals, the SF efficiency v increases with 
the galaxy mass. The IMF is the one of Scalo (1986), and it is 
assumed to be constant in space and time. The use of this IMF is 
motivated by the results by by Kroupa & Weidner (2004), who 
have shown that, since the disc stellar population is made by dis- 
solving open stellar clusters, the disc IMF must be significantly 
steeper than the cluster IMF, because the former results from a 
folding of the latter with the star-cluster mass function. 
An important indication about the IMF in spiral discs comes 
also from chemical evolution models (see Chiappini et al, 1997; 
2001), which clearly indicate that to reproduce the main features 
of the solar neighbourhood and the whole disc of the Milky Way 
a Scalo-like IMF (constant during the entire disc evolution) is 
preferred, and that the Salpeter IMF would overestimate the so- 
lar abundances (see Romano et al. 2005 for a detailed discus- 
sion on this point). In Table 2 we list the adopted values for the 
central surface mass density So, the scaling radius Rq, the star 
formation efficiency v, the infall timescale t,„/ and the IMF for 
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Table 2. Spiral models: parameters 




time/Gyr 



time/Gyr 



Fig. 2. From top left corner, clockwise: SFR, interstellar O abun- 
dance (in units log(0/H)+12), , Type II SNR, and Type la SNR 
vs time for the three spiral galaxy models used in this paper. The 
solid, dotted and dashed lines are the predictions for the models 
with central surface density So = IQQMqI pc^,1Lq - SOOM^jpc^, 
and Eo = IOQQMqI pc^, respectively. 



spirals of different masses. In Fig. |2] we show the time evolu- 
tion of the star formation rates. Type la. Type II SNe rates and 
interstellar O abundance for the three spiral galaxy models used 
in this paper. In spiral discs, to compute the interstellar O abun- 
dance by means of our multi-zone models, we have calculated 
the mass- weighted interstellar O abundance on the whole disc: 



OJSM - 



(5) 



(see Calura & Matteucci 2004), where the o-jsmj are the ISM 
surface densities in the various regions of the disc, whereas the 
Xqj are the gas-phase O abundances in these regions. To com- 
pute Xojs M, for each spiral model, we have considered only the 
regions located at galactocentric distances 3 < R/kpc < 8. This 
is done since the observational metallicities used in this work 
have been obtained considering maximum apertures of ~ 9 kpc, 
as in Savaglio et al. (2005) at redshift z ~ 0.7. On the other 
hand, at the median redshift of 0.07 the SDSS fiber size (3") cor- 
responds to a radius of 3-4 kpc. Our assumption of considering 
only regions at radii < 8kpc has a minor impact on our results. In 
fact, it is worth to stress that, since the innermost regions are the 
densest and most metal enriched ones, the interstellar O abun- 
dance described by eq. |5] is dominated by the contribution of 
regions within R < Akpc. We assume that spiral galaxies do not 
experience galactic winds. This assumption is motivated by the 
results of Tosi et al. (1998), who showed that a chemical evo- 
lution model of the galactic disc including mass outflows is not 
able to reproduce the abundances and abundance ratios observed 
in field stars. 



So 


Rd 


y 




T,„f(R) 


IMF 


{M^lpc^ ) 


(kpc) 


(Gjr-') 




(Gyr) 




200 


2 


0.3 


R- 


1.03 - 1.27 


Scalo 


491 


3.5 


1 


R- 


1.03 - 1.27 


Scalo 


2000 


5 


2 


R 


■0.75-0.5 


Scalo 



Table 3. IiTegular models: parameters 



Baryonic Mass 


Rlum 


V 


Till f 




IMF 


(Mo) 


(kpc) 




(Gyr) 






10« 


1 


0.001 


0.2 


0.25 


Salpeter 


10" 


1 


0.007 


0.2 


0.25 


Salpeter 


lO'o 


1 


0.05 


0.2 


0.25 


Salpeter 



2.3. Irregular galaxies 

For the irregulars we assume a one-zone model with instanta- 
neous and complete mixing of gas inside this zone. IiTegular 
galaxies assemble all their mass by means of a continuous in- 
fall of pristine gas. The SFR is continuous, with SF efficiency 
values lower than the ones used to describe ellipticals and spi- 
rals. In fact, as suggested by Calura & Matteucci (2006a), the 
main features of local galaxies of different morphological types 
may be reproduced by interpreting the Hubble sequence as a se- 
quence of decreasing SF efficiency from early types to late types, 
i.e. from ellipticals to irregulars. 

Also in this case, the star formation can trigger a galactic wind 
if the thermal energy of the gas exceeds its binding energy 
(Bradamante et al. 1998, Recchi et al. 2002). As for Elliptical 
galaxies, for Irr the adopted IMF is the one of Salpeter (1955). 
This choice is in agreement with the results by Calura & 
Matteucci (2004), who have shown that the adoption of a steeper 
IMF in dwarf irregulars leads to an underestimation of their av- 
erage metallicity. 

The rate of gas loss via galactic winds for each element / is as- 
sumed to be proportional to the star formation rate at the time 
t: 



(Gdout = Wi ifr(t) 



(6) 



where w, is a free parameter which describes the efficiency of 
the galactic wind for a given element /, and it is the same for all 
the elements. 

In iiTegular galaxies, the parameter w,-, is fixed in order to repro- 
duce the present-day gas fraction observed in dwarf irregulars. 
The assumption of w, = 0.25 allows us to account for the correct 
fraction of neutral gas in irregulars (Calura & Matteucci 2006a). 
Also in this case, we model irregular galaxies of three different 
luminous masses: M;„,„ = 10**, lO** and lO'^M©. These values 
should bracket the baryonic masses of the most massive typical 
dwarf IiTegular galaxies (Lee et al. 2006), and are compatible 
with the lowest stellar masses of the observational data samples 
used in this work. 

Table 3 shows our adopted values for the initial baryonic mass, 
the luminous radius, the star formation efficiency v, the infall 
timescale and the IMF for each galaxy. For all irregulars, we 
assume a luminous radius of 1 kpc. In Fig. [3] we show the time 
evolution of the star formation rates. Type la. Type II SNe rates 
and interstellar O abundance for the three dwarf galaxy models 
used in this paper. 
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time/Gyr 



time/Gyr 



Fig. 3. From top left corner, clockwise: SFR, O abundance the 
gas (in units log(0/H)+12), Type II SNR, and Type la SNR vs 
time for the three irregular galaxy models used in this paper. The 
solid, dotted and dashed lines are the predictions for the models 
with baryonic mass IO^Mq, lO^M©, and 1O"M0, respectively. 



Case 1) 



t(z,)+A/2 t(z,) 



t(z,)-A/2 



Case 2) 



t(z,)+A/2 t(z,) 



t(z,)-A/2 



Case 3) 

t(zJ 



t(z,)+A/2 t(z,) t(z,)-A/2 



2.4. Different galaxy formation epochs 

As we have just seen, for each galactic morphological type 
we adopt a set of three models correponding to fixed initial 
baryonic masses. These models are designed to reproduce the 
main chemical features of galaxies of different morphological 
types. In general, in chemical evolution studies, the galaxies 
are all considered coeval in the sense that, to compare any 
prediction with local observations, in general one assumes that 
ellipticals, spirals and irregulars have today an age comparable 
to the one of the Universe, i.e. ~ 13.5 Gyr. This is justified by 
the observational fact that the majority of galaxies contains stars 
as old as the Universe. This is certainly true for the ellipticals 
and spirals, although spiral discs are generally younger than 
stellar halos and have ages in the range 8-11 Gyr, whereas for 
the irregulars this is less clear. Then, one considers the model 
outputs at this age, which are compared with observations at 
z = 0. 

However, in a more realistic picture, one should allow for 
some differences among the times at which star formation 
started, in order to take into account a possible intrinsic spread 
in the MZ-relation. For elliptical galaxies, a realistic age forma- 
tion spread is 3 Gyr as indicated by studies of field ellipticals 
(Bernardi et al. 1998), with the cluster ellipticals showing even 
a smaller spread of 2 Gyr (Bower et al. 1992). For the spirals the 
situation is less clear: we tentatively assume a spread of 5 Gyr 
on the basis of the fact that discs take a longer time to assemble 
by gas accretion than spheroids. As we will see later (Sec. l3.4l l. 
such an age dispersion values are compatible with other age es- 
timates for spirals present in the literature. 
Concerning irregulars, in principle, the age spread could be as 
long as ~ 10 Gyr, i.e. comparable to a Hubble time. 
In this section, we describe how we simulate a population of 
non-coeval galaxies and how we create a finer mass grid. First 
of all, we define as galaxy formation epoch the time at which 
star formation starts in a given galaxy. For each morphological 
type, we choose a mean redshift of formation z/ and an age dis- 
persion A,, relative to the spread in the galaxy formation time. 
If the mean redshift of formation is Zf and t{z) is the age of the 
Universe at a given redshift z, galaxies are allowed to form across 
the cosmic time interval t{zf) - l^tl'^ < t < t(zf) + At/2, i.e. this 
means that galaxies may start forming at a redshift greater than 
Zf. If Zo is the redshift of observation of a particular galactic 
observable, the ages (defined as the times elapsed since the be- 
ginning of star formation) of the galaxies at Zo are allowed to 
vary between two values r„„„ and T,„a_^, depending on Zo and the 
assumed values of Zf and A,: 

- If Hzf) - A,/2 < t(zo) < t(zf) + At/2, i.e. if galaxies are still 
forming at the redshift of observation, r„„„ - and T„t„^ = 
abs(t(zf) - (t(Zo) + At/2) (Cases 1 and 2 in Figgli 

- If t{zo) > t(zf) + At/2, i.e. galaxies have stopped forming at 
the redshift Zo, then = t(zo) - (f(z/) + At/2) and T,„ax = 
t(Zo) - (t(zf) - A,/2) (Case 3 in Fig© 



Fig. 4. In cases 1 and 2, galaxies are still forming at the redshift 
of observation, i.e. t{zf) - A, < t{zo) < t{zf) + A,, where f(z) is 
the age of the Universe at the redshift z- In case 3, galaxies have 
stopped forming at the redshift Zo, so t{zo) > t{zf) + A,. 



where t(zf) and t(zo) are the ages of the Universe at the forma- 
tion redshift Zf and at the observation redshift Zo, respectively. 
For each morphological type and for any given mass, we ex- 
tract random ages from a flat probability distribution between 
Tmi„ and T,„„^ and we compute the relevant physical quantities at 
this randomly extracted age. With our method, we can simulate 
a continuous galaxy formation process in a given time interval 
and a finite dispersion of the galaxy formation epoch. At any 
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redshift, this method is useful to generate an almost-continuous 
galactic stellar mass grid. 

In Fig.|5] we show an example of the application of this tech- 
nique: the predicted MZ relation for ellipticals, spirals and irreg- 
ulars, obtained assuming zy = 3, A, = 3 Gyr and computed at a 
redshft Zo - 2.2. For each morphological type, the black curves 
represent the evolutionary tracks of the baseline models of three 
baryonic masses, whereas the grey regions represent the MZ re- 
lation computed by extracting randomly stellar masses and O 
abundances at times between r,„,„ = Gyr and T,„ax - 2.3 Gyr. 
The discontinuity predicted for elliptical galaxies (lower panel 
of Fig. |5]l is due to the steep increase of the metallicity in the 
first evolutionary phases, i.e. within the first 0.1 Gyr (see Fig.[T]). 
Also the evolution of the stellar mass density is fast, hence small 
variation in the extracted times correspond to large variations for 
both the mass and the metallicity. 

The method described in this section will be used to gener- 
ate mass-metallicity plots for galaxies of various morphological 
types, to be compared with observational data obtained at vari- 
ous redshifts. 

3. Results 

In Fig. |6] we show the predicted MZ-relations at the present 
time for spirals and irregulars and compare them with the best 
fit of the MZ-relation obtained by Maiolino et al. (2008) from an 
analithical fit to the local observations (Kewley & Ellison 2008), 
and to the data of Lee et al. (2006), who determined the MZ 
relation in a sample of local dwarf galaxies. The agreement be- 
tween our predictions and the observations is quite good. It is 
worth to stress that we do not show the results for local ellip- 
ticals since in this paper we are comparing gas abundances at 
various redshifts in star forming galaxies, so we cannot compare 
the predicted present time MZ-relation for the gas in ellipticals 
with the data, since these galaxies stopped forming stars several 
Gyrs ago. On the other hand, local ellipticals show a definite 
MZ-relation, based on the abundances of their stellar popula- 
tions, which is well fitted by our models (see Pipino & Matteucci 
2004). 

Figure |6] shows that the local MZ relation can be naturally 
explained by assuming a lower star formation efficiency for 
less massive galaxies, irrespective of the galactic morphologi- 
cal type, and that galactic winds are not necessary to explain it. 
In fact, our models for spirals do not include galactic winds. The 
reason for this choice is that galactic winds are not required to 
explain the main features of disc galaxies, as shown by Tosi et 
al. (1998). We consider galactic winds in dwarf irregulars, since 
they are observed in these systems. However, these winds do not 
carry away large amounts of matter, in agreement with the fact 
that irregulars are gas-rich systems. 

It is worth to note that Dalcanton (2007) reinterpreted the ob- 
served MZ-relation in terms of infall/outflow in galaxies by 
means of simple models with IRA and concluded that only gas 
rich systems with low star formation rates can produce and main- 
tain low effective yields. This result shows the importance of as- 
suming a lower star formation efficiency in lower mass systems. 
The flattening of the MZ relation observed at M» = IO^'Mq is 
accounted for by our models and is due to the assumed SF effi- 
ciencies of the high mass spiral disc models. The adopted SF effi- 
ciencies cause a steep evolution of the (O/H) vs time relation (see 
Fig.|2]). Note that for larger galaxies, their present-day metallic- 
ity is reached at earlier times . In fact, for the MlOl model, very 
little evolution in the (O/H)- time diagram is predicted at cos- 
mic times greater than ~ 5 Gyr The growth of the metallicity 



for the MW model is slightly slower than the one predicted for 
MlOl. This is due to the fact that the adopted SF efficiencies are 
comparable. On the other hand, the model for small discs has 
a SF efficiency lower by almost one order of magnitude, hence 
the growth of metallicity with time is slower, and the present- 
day metallicity is much smaller than the ones of the MW and of 
MlOl. 

In conclusion, the observed flattening of the MZ relation for 
large discs indicates that in the past, their SF efficiences must 
have been similar This is consistent also with the fact that the 
O gradients observed in the MW and in MlOl are very similar, 
although MlOl is more massive (see Chiappini et al. 2003). 

3. 1 . The calibration of tiie mass-metallicity at redshift zero 

In Fig. |7] we show the predicted MZ relation for local star 
forming galaxies, compared with a set of observational MZ 
relations at z = 0.07. In our models, only spirals and irregulars 
experience star formation at z = 0.07. This is supported by 
several observational studies, which have shown that local 
ellipticals are dominated by old stellar populations (de Freitas 
Pacheco, Michard, & Mohayaee 2003 and references therein). 
The filled contours in Fig. |7] are the predicted MZ relation 
computed at z=0.07 for spirals and irregulars, whereas the 
lines are fits to the observed local mass-metallicity relation 
as reported by Kewley & Ellison (2008), obtained assuming 
different metallicity calibrations. 

Kewley & Ellison considered ~ 28000 star-forming galaxies of 
the SDSS with derived stellar masses. They applied 10 different 
metallicity calibrations to the emission lines of the spectra of 
this sample, in order to investigate the effect of the calibration 
on the MZ relation. Some calibrations are based on a direct 
method for determining metallicity. This method, often referred 
to as the 'T(." method, is based on the measure of the ratio of 
two lines, such as the 04363 auroral line and the [OIII]/l 5007 
lines and allows an estimate of the electron temperature of the 
interstellar gas, finally used to determine the metallicity. Some 
other calibrations are empirical, and are obtained by fitting 
the relation between direct metallicities and strong line 
ratios for Hji regions, such as [NII]/16584/Hq'. Finally, in other 
cases, mostly at high metallicity, where Tg is not measurable, 
photoionization models are used to compute strong-line ratios. 
Is it important to stress that each of these calibrations are valid 
in relatively narrow metallicity intervals (Kewley & Ellison 
2008). 

In Fig. U\ the curves labeled KK04 (Kobulnicky & Kewley 
2004), KD02 (Kewley & Dopita 2002) are based on theoretical 
methods, the curves labeled PP04 (Pettini & Pagel 2004) and 
P05 (Pilyugin & Thuan 2005)) are based on empirical methods, 
while the D02 (Denicolo, Terlevich & Terlevich 2002) curve 
has been derived by means of a combined method. The sample 
of D02 is composed by a set of Hji regions, some of which have 
metallicities derived using the method, while some others 
have empirical or theoretical metallicities. 
As can be seen in Fig.|7] the choice of the metallicity calibration 
plays an important role in determining observationally the MZ 
relation. Kewley & Ellison (2008) have outlined that different 
calibrations can produce very different results, concerning both 
the zero-point and the slope of the observational MZ relation, 
as well as the location of the plateau for high-mass galaxies. 
For the particular cases reported in Fig. [T] including also the 
most extreme calibrations, one can see that the zero point of 
the observed MZ relation can vary by ~ 0.4 dex, whereas the 
high-mass plateau position may vary by ~ 0.3 dex, without 
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Fig. 5. Predicted MZ relation for galaxies of various morphological types. In each panel, the curves represent the evolutionary tracks 
of individual galaxies of various baryonic masses. The grey areas are the predictions obtained by means of our method to create an 
almost-continuous galaxy formation and a finer stellar mass grid (see Sect. 12.4b . In this case, we have assumed z/ = 3, A, = 3 Gyr 
for all galaxies. The predicted MZ is computed at redshft Zo = 2.2:This choice is due to the fact that at this cosmic epoch also the 
ellipticals are still star forming galaxies. 



considering the extreme calibration by Pilyugin & Thuan 
(2005), leading to an almost flat MZ. 

Concerning the predictions, here we have assumed z/ = 3 
and A, = 5 Gyr for spirals, and zy = 3 and A, = 10 Gyr 
for irregulars. However, as we will see later, for spirals and 
irregulars the parameters z / and A, have little effect on the zero 
point and on the slope of the predicted mass-metallicity relation, 
whereas they mostly influence the predicted dispersion of the 
stellar masses and of the O abundances. As discussed in Sect. |2l 
the models used for spirals and irregulars are calibrated in order 



to reproduce a large set of chemical evolution constraints from 
local observations, such as the abundances observed in stars of 
the Milky Way, the stellar metallicity distribution in the solar 
neighbourhood (Chiappini et al. 2001), the abundance gradients 
in the MW and in MlOl (Chiappini et al. 2003), as weU as 
the stellar abundances and stellar metallicity distributions in 
local dwarf galaxies (Lanfranchi & Matteucci 2003). Since our 
models can reproduce such a large amount of local independent 
observations, it seems appropriate to use them to find some 
constraints on the best calibration method for the local observa- 
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tional MZ relation. Among the different calibrations considered 
here, our predictions are compatible with the MZ relations 
obtained by adopting the calibrations KD02, D02 and PP04. 
The possibility to use our models to constrain the metallicity 
calibrations is motivated by the fact that most of the observables 
cited above and used to tune our models have typical errors 
much lower than the calibration uncertainty. For instance, 
the abundances observed in local stars have typical errors 
of 0.05-0.1 dex, in extremely rare cases larger than 0.2 dex 
(Cayrel et al. 2004; Spite et al. 2005). Similar uncertainties are 
typical also for abundances observed in the stars of local dwarf 
spheroidal galaxies (Shetrone et al. 2003) and blue compact 
galaxies (Izotov & Thuan 1999) and in local early type galaxies 
(Thomas et al. 2005). 

Concerning the calibration method used by Lee et al. (2006) for 
the data reported in Fig.|6] for most of the objects the method 
was used, consistent with the one of Maiolino et al. (2008) for 
the metallicity range investigated by Lee et al. (2006). For a very 
limited number of systems, Lee et al. (2006) used the strong line 
method. Given the uncertainties of various calibration methods 
in particular at metallicities logiOjH) +12 < 8.2, we neglect 
this difference and use the data by Lee et al. (2006) at face value. 



3.2. The observational data used in ttiis work 



The observational data for the MZ relation at redshits z=0.07, 
z=0.7, z=2.2 and z=3.5 are taken from Kewley & Ellison (2008), 
SavagUo et al. (2005), Erb et al. (2006a) and Maiolino et al. 
(2008), respectively. These are all MZ relations for the gas phase 
metallicity, and include only star forming galaxies (active star 
formation is required to ionize the gas which produce the nebular 
lines used for the metallicity determination). In these works dif- 
ferent diagnostics and calibrations have been adopted, which af- 
fect the zero point and slope of the MZ relation. This issue is dis- 
cussed extensively by Kewley & Ellison (2008), who also pro- 
pose a criterion to interlace different calibrations. However, each 
of these calibration is appropriate in different, relatively narrow 
metallicity ranges. Maiolino et al. (2008) note that the metallicy 
range spanned by galaxies through the cosmic epochs is so wide 
(7.7<12+loh(0/H)<9.1) that no single metallicity calibration 
method is appropriate to cover it. As a consequence, Maiolino 
et al. (2008) define a new calibration that adopts the electron 
temperature method at low metallicities (12+log(0/H)<8.3) and 
the photionization models by Kewley & Dopita (2002) at high 
metallicities (12-i-log(0/H)<8.3). Details on this method and the 
appropriateness of this calibration scale are given in Maiolino et 
al. (2008). To have the metallicities of all galaxies at the various 
redshifts sampled by the different studies on a consistent metal- 
licity scale, Maiolino et al. (2008) re-determined the metallici- 
ties of the previous works at lower redshifts by using the new 
calibrations inferred by them. Maiolino et al. (2008) also correct 
the stellar masses of different studies to the same IMF. In the 
following section, we will use the same data as the ones used 
by Maiolino et al. (2008), together with SFR determinations at 
various redshifts. As we will see later, the combined use of ob- 
servational mass-metallicity measured and of SFR determina- 
tions will allow us to put some constraints on the morphologies 
of the galaxies populating the observed MZ relation at various 
redshifts. 



3.3. Tlie evolution of tlie mass-Metailicity relation 

In Figure [8] we show the observed redshift evolution of the MZ 
relation in star forming galaxies, as well as the predicted MZ 
relation for various galactic morphological types used in this 
work. The predicted MZ relations for ellipticals, spirals and ir- 
regulars have been computed by assuming zf - 3, A, =3 Gyr 
The black points are the observations at various redshifts, de- 
scribed in Sect. 13.21 whereas the blue and red contours represent 
the regions where 100% and 75 % of our predictions lie, respec- 
tively. 

From this plot, it is clear how galaxies of various morpho- 
logical types occupy different regions of the MZ plot. At any 
redshift, in general, elliptical galaxies present the highest stellar 
masses and the highest metallicities, whereas the irregulars are 
the least massive galaxies, characterized by the lowest O abun- 
dances. In this case, ellipticals appear only at z - 3.5 and at 
z = 2.2 since, with the age dispersion of 3 Gyr as chosen here, 
these galaxies become passive at redshift < 1.4. From Fig. [8] we 
can see that the majority of ellipticals show metallicities higher 
than the observations, both at z = 3.5 and at z = 2.2. This 
fact concerns only ellipticals and may be due to the fact that the 
growth of (O/H) vs time is very steep for elliptical galaxies. We 
will investigate this issue in deeper detail later, in Sect. 13.51 On 
the other hand, spiral galaxies show metallicities always compa- 
rable to the values observed at various redshifts. Irregular galax- 
ies have low star formation rates and at high redshift their stel- 
lar masses are very small, in general M, < lO^M© at z > 2.2. 
The observed masses of the galaxies populating the MZ rela- 
tion at these redshifts are in general higher than M, ~ lO^M©, 
hence there is poor overlap between the observational MZ rela- 
tions and our predictions for irregulars at z > 2.2. At z = 0.07, 
the SDSS catalogue contains also a lot of galaxies with stellar 
masses M, < IO'Mq. The galaxies with these stellar mass val- 
ues overlap well with our predictions for irregular galaxies. At 
z < 0.7, we note also that the chosen value for the age disper- 
sion does not allow us to reproduce the observed data dispersion. 
Also this aspect will be investigated in deeper detail later on. 
A comparison between Fig.[8]and Fig.|9]is useful to appreciate 
the effects of our parameter Z/ on our predictions for the MZ 
relation as a function of redshift. The results for spirals and ir- 
regulars are fairly insensitive on the assumption of the redshift 
of formation z /, whereas for ellipticals the assumption of an high 
formation redshift, such as z/ = 6 in Fig.|9] leads to a population 
of extremely high-metallicity ellipticals at z = 2.2, which does 
not match the MZ relation observed at this redshift by Erb et al. 
(2006a). 

It is maybe worth stressing that the contours are not linked to the 
comoving densities of galaxies of various morphological types 
present at any redshift. Such a prediction would require the adop- 
tion of a galaxy luminosity function or a stellar mass function as 
a function of redshift. In this case, the predictions represent the 
number of galaxies present in each M-Z bin, normalized to the 
total amount of simulated galaxies. The blue contours enclose 
the region where at least I galaxy is present. On the other hand, 
the red contours enclose the region where the total number of 
galaxies is 0.75 times the total number of simulated galaxies. 
In order to relie on a robust statistics, we have simulated ~ 10^ 
galaxies, i.e. performed ~ 10^ interpolations on the mass and 
metallicity grids according to the method described in Sect. 12.41 
However, we have verified that in most cases, the shape of the 
contour plots is not strongly dependent on the adopted number 
of simulated galaxies. The adopted number of simulated galax- 
ies may have some effects in regions of the plots where galaxies 
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are rare, correspoding to the discontinuities visible in the con- 
tours. An example is the peculiar "comb" feature present in the 
M-Z plots of elliptical galaxies in Fig. [8] on the left bottom side 
of the plot. 

As shown by Calura et al. (2009), the use of a cosmological 
galaxy formation model can provide directly galaxy abundances 
in the MZ plot, but no indication concerning the galaxy mor- 
phology. In the following, we will present a method to constrain 
the morphology of galaxies, starting from the assumption that 
morphology does not change with redshift. 

3.4. The mass-metallicity relation of star-forming galaxies: 
constraints on their morphology 

In this section, we focus on the observed evolution of the mass- 
metallicity relation for star-forming galaxies, considering also 
their star formation rates. The SFR provides us with a funda- 
mental information, since, as we will see, it allows us to have 
further hints on the nature of the galaxies building the MZ rela- 
tion. Furthermore, a study of a 3-dimensional plot linking mass, 
metallicity and SFR is a crucial test to our models, which should 
be able to reproduce at the same time all of these properties at 
any redshift. 

In FiglTO] we show the observed and predicted MZ, (0/H)-SFR 
and SFR-Mass plots for all morphological types as a function of 
reshift. At each redshift, the observed quantities are compared 
to theoretical predictions, obtained for ellipticals, spirals and ir- 
regular galaxies as explained in sect. 12.41 For the observational 
SFRs, at z = 0.07 we use the values derived for SDSS galaxies 
by Brinchmann et al. (2004). At z = 0.7 we use observations 
by Maier et al. (2005) and Juneau et al. (2005), at z = 2.2 we 
use the values by Erb et al. (2006b) and at z = 3.5 we use the 
data by Maiolino et al. (2008). It is important to note that, at 
z - 0.07, z = 0.7 and z = 3.5 we use metallicities, stellar masses 
and SFRs estimated observationally for individual galaxies. In 
particular, at z = 0.7, we use metallicity measurements only for 
those galaxies whose SFRs are available. At z = 2.2, the mea- 
surements of O/H for the individual galaxies are not available. 
The only data available are represented by the rebinned MZ re- 
lation as published by Erb et al. (2006a). On the other hand, 
the SFRs for individual galaxies are available (Erb et al. 2006b). 
For these galaxies, we calculate the metallicity from the stellar 
mass, by means of the relation found by Maiolino et al. (2008) 
for star-forming galaxies at z = 2.2, which basically represents 
an analithical fit to the data of Erb et al. (2006a). 
For the parameters z / and A,, we have investigated several cases, 
in the ranges 2 < z/ < 6 and 2 Gyr < A, < 10 Gyr. In Fig. [TOl 
we show the case computed with the fiducial set of the parame- 
ters, i.e. Z/ = 3 for all galaxies. A, = 3 Gyr for elliticals. A, = 5 
Gyr for spirals and A, - 10 Gyr for irregulars. It is worth to 
note that for irregulars and spirals, the assumption of a smaller 
A, does not change the shape of the MZ relation, of the O/H 
vs SFR and of the SFR vs M, relations, but affects only the pre- 
dicted dispersion. This plot shows how most of the observational 
constraints considered in this work are satisfactorily reproduced 
once all morphological types are included. In this figure, the pre- 
dictions for ellipticals, spirals and irregulars are shown with dif- 
ferent colours. At z = 3.5, the observed MZ, metallicity vs SFR 
and SFR vs stellar mass plots are well reproduced by our models 
for elliptical galaxies, whereas the spirals and irregulars present 
SFRs considerably lower than the observations. This indicates 
that the galaxies observed by Maiolino et al. (2008) are likely 
proto-ellipticals observed during the starburst phase. 
At z = 2.2, a morphological mix of spirals and ellipticals can 



reproduce the MZ relation, the (0/H)-SFR and the SFR-Mass 
plots as observed by Erb et al. (2006a, 2006b). However, the 
figure indicates that there is a limited number of observational 
points which are compatible with both predictions for ellipticals 
and spirals. This leads us to a degeneracy concerning the inter- 
pretation of the observational results. A possible way to break 
this degeneracy will be suggested later in Sect. 13.61 
At z = 0.7, by means of the models for spirals we slightly under- 
estimate the metallicities in the MZ and in the (O/H) vs SFR 
plots, however, given the uncertainties due to calibration dis- 
cussed in Sect. 13.11 we do not consider this as a major issue. 
On the other hand, in the SFR vs mass plot, we predict a posi- 
tive correlation between M» and the SFR. The observational data 
seem to suggest a weaker correlation between M, and the SFR. 
Furthermore, several observed SFRs are higher than the values 
we predict for spiral galaxies. This discrepancy could be partly 
reduced by considering larger age dispersions for spirals. For 
spirals, assuming a redshift of formation z/ = 3, an age disper- 
sion of A, = 5 Gyr imphes present ages between 8.86 Gyr and 
one Hubble time. Independent estimates of the present ages of 
spiral galaxies indicate values ranging from 5-6 Gyr up to one 
Hubble time (Bell & de Jong 2000; Boissier et al. 2001), with 
low mass spirals younger than high mass discs. This indicates 
that the assumed age dispersion for spirals may be slightly un- 
derestimated, in particular for low mass spirals. However, the ob- 
served SFRs may be higher owing to episodic starbursts, possi- 
bly triggered by a dynamical process such as galaxy interactions 
(Alonso-Herrero, Rieke & Rieke 1998). These events, which are 
of random nature and which are not taken into account by our 
models, could be the cause of the weak correlation between SFR 
and M, as observed at z ~ 0.7. For this reason, this explanation 
may be the most likely for the discrepancy between our predic- 
tions and the observed SFRs at z = 0.7. 

At z = 0.07, we can reproduce the shape of the observed MZ 
relation. However, the dispersion of the data is larger than what 
our predictions indicate. For spirals, by assuming A, = 5 Gyr, 
we predict a maximum dispersion of 0.2 dex. The observed 
dispersion cannot be reproduced even assuming the unrealis- 
tic value of A, = 14 Gyr, i.e. equal to an Hubble time. Also 
in this case, this discrepancy could be due to the fact that the 
star formation histories considered here for spirals and irregu- 
lars do not take into account stochastic events such as episodic 
starburtsts, which could increase the dispersion of the observed 
metallicities. However, this does not represent a major concern 
since, as shown in Kewley & Ellison (2008), also the dispersion 
depends on the calibration. The use of the metallicity calibra- 
tion of KD02, along with some others (e.g. Pilyugin & Thuan 
2005) produce dispersions in the local MZ relation of ~ 0.6 dex, 
whereas with other calibrations the degree of dispersion in the 
data is lower (~ 0.4 dex or lower,Kewey & Ellison 2008). For 
this reason, in this study, the observed y-axis dispersion in the 
MZ or (O/H) vs SFRs cannot be considered as a key-constraint 
to our models. 

Also at z = 0.07, a large number of observed galaxies present 
SFRs higher than what our predictions for spirals indicate. 
However, the observed correlation between SFR and M, is here 
reproduced by our predictions. Interestingly, at z = 0.07 the 
data indicate SFRs even higher than those observed at z = 0.7. 
Locally, very high SFR may be measured in luminous infrared 
galaxies or in galaxies experiencing episodic starbursts, not con- 
sidered here. However, it is possible that, for some galaxies, the 
star formation rates may be overestimated. In fact, by compar- 
ing to other data of star-forming galaxies observed in the AEGIS 
field (Noeske et al. 2007), in the lowest redshift bin the max- 
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imum observed SFRs are log{S FR/MQ/yr) ~ 1.3, consistent 
with our values predicted for spirals, whereas the SFRs of the 
SDSS data may be as high as log(S FR/Me/yr) ~ 1.8 or, in a 
very few cases, even higher. We suggest here that the discrep- 
ancy may be due to different dust extinction corrections. This as- 
pect will be investigated in detail with next release of the SDSS 
data. 

3.5. A possible dust-obscuration bias atz- 3.5 

In Sect. 13.31 when discussing the MZ relation at z = 3.5, we 
have seen that our predictions for ellipticals, i.e. the best can- 
didates for the galaxies observed at z = 3.5 by Maiolino et al. 
(2008), indicate metallicity values higher than the observed ones 
by ~ 0.5 dex. This difference is higher than the offset due to 
different calibrations discussed in Sect. 13.11 since the metallicity 
calibration we use in this paper provides metallicities very close 
to the ones of Kewley & Dopita (2002), one of the calibrations 
providing the highest metallicities in the set analyzed by Kewley 
& Ellison (2008). This fact concerns only the predictions for el- 
lipticals galaxies. The discrepancy between the metallicities pre- 
dicted for most ellipticals at z = 3.5 and the observed values may 
be partly due to dust obscuration. To investigate this issue, we 
use chemical evolution models for ellipticals including also dust 
grain production and destruction. The method to model dust evo- 
lution in galaxies is described in Calura, Pipino and Matteucci 
(2008a). Dust grains are produced in low mass stars. Type la 
and Type II SNe, destroyed by SN shocks and the grain cores 
are allowed to accrete mantles, i.e. to grow by mass, during the 
starburst, when large reservoirs of molecular H is available for 
star formation. For further details on the formalism built to han- 
dle dust evolution in galaxies of different morphological types, 
we refer the reader to the paper by Calura et al. (2008a). 

In this case, we are interested in the dust surface density at 
various evolutionary phases. In ellipticals, the dust surface mass 
density cr^ is 

CTd = MdlnRlff (7) 

where is the total dust mass, including the refractory 
elements C, O, Fe, S, Si, Fe and Mg (Calura et al. 2008a), and 
Reff is the effective radius. 

A study of the amount of dust present in proto-ellipticals at 
various evolutionary stages is interesting because of the link be- 
tween dust mass and obscuration (Calzetti 2001). Stages char- 
acterized by higher dust surface densities coiTespond to more 
dust-obscuration, hence in principle to phases more difficult to 
observe. However, as stressed before, owing to the intense UV 
radiation field present in //// regions, dust depletion is not likely 
to affect the measured metallicities (Okada et al. 2008). 
In Fig.[TT] we show the predictions for the MZ plot and the SFR 
vs Mass plot for ellipticals at z = 3.5 grouped in three differ- 
ent regions, depending on the dust column density at various 
evolutionary stages. The data points with error bars are the ob- 
servations of Maiolino et al. (2008). The three different colours 
of the models indicate predicted galaxies with cr^ < O.IM©/ pc- 
(magenta), 0.1 < a-dlM^lpc^ < 0.3 (cyan), crj > 0.3Mq/pc^ 
(green). The solid lines represent the single evolutionary tracks 
for the three models used to describe elliptical galaxies. The 
dashed lines are isochrones at times 0.05 Gyr, 0.1 Gyr and 0.2 
Gyr. 

In the MZ plot, most of the data points overlap with the predic- 
tions computed for ellipticals at the earlier phases, characterized 



by the lowest dust surface densities. The fact that some stellar 
masses are underestimated by our predictions may be due to the 
uncertainties affecting the spectro-photometric models used to 
determine the observational M,. In this case, we plot the stellar 
masses obtained by Maiolino et al. (2008) by adopting the spec- 
tral synthesis models of Bruzual & Chariot (2003). Maiolino 
et al. (2008) published also stellar masses obtained with the 
Maraston (2005) spectral synthesis models, in general produc- 
ing stellar masses lower by ~ 0.1 dex or less. The use of these 
stellar masses would reduce the discrepancy between our predic- 
tions and the observations. 

In the SFR-Mass plot, most of the observations, plotted with 
their error bars, overlap with predictions for galaxies with inter- 
mediate dust surface densities. An important caveat concerning 
the observed SFRs we should keep in mind is that, in most cases, 
it is likely that the observed SFRs may represent underestima- 
tions to the real values. In general, in observational studies of 
high-redshift galaxies, the SFRs and the ages are very uncertain 
parameters. The SFRs are determined from the UV luminosity 
by assuming that the age of each galaxy is larger than 50 Myr. 
By removing this constraint, the estimated SFRs may be higher 
by a factor of 3 for the youngest objects. As a consequence, all 
the points would move upwards, towards regions characterized 
by lower dust surface density, hence less dust-obscured. 
The evolution of the ellipticals is very quick within the first 0.2 
Gyr (see fig. [TJ, such as the accumulation of a large dust mass. 
Fig, nn indicates that, if the AMAZE sample at z=3.5 is com- 
posed primarily by the progenitors of ellipticals, most of these 
galaxies may remain observable for times < 0.2 Gyr, consider- 
ing also the uncertainties in the observed SFRs. After 0.2 Gyr, 
the bulk of galaxies may become heavily dust-obscured and can- 
not be detected by current surveys in the optical-UV bands. 
A population of extremely dust obscured and vigorously 
star-forming galaxies is represented by the SCUBA galaxies 
(Clements et al. 2008), very luminous in the far infrared and 
sub-millimetric bands owing to massive dust reprocessing, but 
intrinsically faint in the optical band. Schurer et al. (2009) have 
used the same elliptical models described in this work to com- 
pute the photometric properties of galaxies, by taking into ac- 
count both the chemical and the spectro-photometric evolution 
of silicate and carbon dust grains. Schurer et al. (2009) showed 
that the general shape of the spectral energy distribution and the 
observed amount of dust reprocessing of SCUBA galaxies are 
correctly reproduced by means of our models. 
In summary, for ellipticals we predict metallicities larger than 
the observation owing to the fact that the evolution of our galax- 
ies is very fast. A possible way to reduce the predicted metallic- 
ities could invoke longer infall times for the ellipticals. A longer 
infall timescale would cause in each galaxy a longer star for- 
mation period and lower star formation rates, consequently a 
slower growth of the metallicity versus time. Another possibility 
could invoke stronger outflows in the elliptical models, or metal- 
enhanced ouflows, where the metallicity of the ejected matter is 
higher than the one of the ISM. However, this would require the 
introduction of additional parameters in the elUptical model. 

3.6. MZ relation for various elements and interstellar 
abundance ratios 

In Fig. [121 we show our predictions for the MZ relations for 
spiral galaxies, considering on the y axis various chemical el- 
ements: O, C, N, Mg, Si, Fe. Here we have assumed Zf - 3 
and A, = 5 Gyr. These elements are produced by different stel- 
lar sources and restored into the ISM on different timescales 
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(see Matteucci 2001; Calura & Matteucci 2006b). C and N are 
produced mainly by low and intermediate mass stars of mass 
< 8Mq. Mg, Si and O are produced by massive stars, dying as 
Type II SNe. However, Mg and Si are produced also by Type la 
SNe in a non negligible amount. Finally, the bulk of Fe is pro- 
duced by Type la SNe. The maximum dispersion of the predicted 
MZ relations (i.e. the dispersion at stellar masses lO^'^ - lO^'^M© 
) reflect the different production timescales of the elements and 
the lifetimes of the stars where they are synthesised. The small- 
est dispersion is of ~ 0.2 dex, visible for O. Dispersions of 
~ 0.22 - 0.25 dex are predicted for C, N and Mg. The largest 
dispersions are predicted for the elements Si (~ 0.3 dex) and Fe 
(~ 0.4 dex), produced in considerable amounts by Type la SNe. 
This is related to the fact that the lifetimes of stars exploding 
as Type la SNe span a very large range, from 0.03 Gyr for the 
highest masses (M = 8Mq ) up to a Hubble time for the low- 
est mass stars. The reason for the larger scatter in the elements 
produced by Type la SNe with respect to the a elements can be 
understood by looking at Fig.|2] While the Type II SN rate, scal- 
ing with the rate of production of the a elements, presents a peak 
at early times and then decreases progressively, the Type la SN 
rate presents a monotonic increase with cosmic time, hence the 
massive production of the Fe-peak elements lasts for a longer 
time than the a-element production. 

3.7. The importance of abundance ratios 

In Sect. |33] we have seen that, in the MZ, O/H vs SFR and 
SFR vs Mass plots, the predictions for ellipticals and spirals 
overlap, thus making diflicult the interpretation of the obser- 
vational data. A way to overcome this issue is to consider di- 
agrams which involve abundance ratios between two different 
chemical elements. The abundance ratios between two chemi- 
cal elements synthesized on different timescales can be used as 
"cosmic clocks", providing us with information on the relative 
roles of various stellar sources in the chemical enrichment of the 
interstellar medium (Matteucci 2001). In particular, the study of 
the [ff/Feflis of major importance, owing to the difference in the 
timescales for a-elements and Fe production. 
In Fig. [13] we show the predicted [O/Fe] plots computed at 
Z - 3.5 and at z = 2.2 for ellipticals, spirals and irregulars. 
As one can see, in both plots the three morphological types oc- 
cupy very distinct regions. This eff'ect is primarily due to the star 
formation histories: ellipticals galaxies form quickly, by means 
of intense starbursts. For this reason, at low metallicity (e.g. 
[Fe/H]< -1) the enrichment of their ISM is dominated by a- 
elements, produced by short-living massive stars. Owing to their 
short and strong starbursts, ellipticals reach high metallicities 
very quickly by means of Type II SNe, thus having high [a/Fe] 
ratios for a large range of [Fe/H] . On the other hand, spirals and 
irregulars have lower star formation rates and more prolonged 
star formation histories. For this reason, they reach high metal- 
licities at later times, when the contribution of Type la SNe, pro- 
ducing the bulk of Fe and exploding on long timescales, becomes 
important. In principle, if measures of interstellar Fe abundances 
in high-redshift star-forming galaxies were accessible, thanks to 
our chemical evolution models, it would be possible to have im- 
portant hints on the morphology of these objects. 



' All the abundances between two different elements X and Y are 
expressed as [X/Y] = log(X/Y) - log(X/Y)o, where (X/Y) and (X/Y)o 
are the ratios between the mass fractions of X and Y in the ISM and in 
the sun, respectively. We use the set of solar abundances as determined 
by Grevesse at al. (2007). 



The measure of metal abundances in high-redshift star forming 
galaxies has become accessible in the recent years, thanks to 
the use of 10 m telescopes for spectroscopic studies, for which 
bright Lyman-break galaxies observed at z > 2 are particu- 
larly valuable tools. Unfortunately, in these objects the measure 
of abundances other than O are hampered by weak emission 
lines, depletion uncertainties and unknown ionization correc- 
tions (Pettini 2004; Leitherer 2005). More promising objects for 
interstellar abundance ratio studies are Damped-Lyman Alpha 
galaxies, observed in absorption in the spectra of background 
Quasars. For these systems, the use of high-resolution spectro- 
graphs offers the access to weak interstellar lines, which lead to 
precise determination of interstellar abundances for various ele- 
ments (Dessauges-Zavadsky et al. 2007). Recently, an attempt to 
study the possible existence of a mass-metallicity relation also 
in Damped Lyman Alpha systems has been suggested in a pi- 
oneeristic work by Ledoux et al. (2006). The interpretation of 
these data by means of models, such as the ones presented in 
this paper, will represent the next step in cosmic chemical evo- 
lution studies. 

4. Conclusions 

In this paper, we have performed a study of the evolution of the 
MZ relation in galaxies of different morphological types, and 
compared our results with the most recent observational data. 
In particular, we have used various observations achieved by 
several different authors, adopting a common metallicity cali- 
bration suited to our study, where the measured metallicities of 
the galaxies, observed at redshifts from 0.07 to 3.5, span a large 
range of values. We have used chemical evolution models for 
ellipticals, spirals and irregular galaxies. The models used here 
are able to reproduce the main properties of local galaxies. For 
each morphological type, we have used three baseline models 
of different present-day stellar masses. We have assumed that 
galaxy morphologies do not change with cosmic time. We have 
developed a method to take into account a possible spread in the 
epochs of galaxy formation, namely in the times at which differ- 
ent galaxies start forming stars. 

Our main results can be summarized as follows. 
1) Previous papers interpreted the observed MZ-relation in 
terms of the Simple Model of galactic chemical evolution. 
In this model, the main quantity linked to the galactic metal 
content is the "effective yield" yz - ZjlniiT^), where yu is the 
gas fraction. In reality, galaxies will suffer infall/outflow and 
the true yield will be lower than the effective one. The strongest 
assumption about the Simple Model is the instantaneous recy- 
cling approximation, and if elements different than O would 
be considered, such as N or Fe, this hypothesis would be very 
poor, since both N and Fe are produced mainly by low and 
intermediate mass stars on long timescales. For O, the IRA can 
still be acceptable. Therefore, interpreting the MZ-relation in 
terms of the Simple Model has led several authors to conclude 
that larger galaxies should have a larger true yield than less mas- 
sive ones. This can be achieved in several ways: i) by varying 
the IMF, ii) by decreasing the importance of outflows in more 
massive galaxies and iii) by decreasing the importance of infall 
in more massive galaxies. In this paper, we have shown that a 
MZ-relation naturally arises if less massive galaxies experience 
a lower star formation efficiency (i.e. the star formation rate per 
unit mass of gas), irrespective of the galactic morphological 
type. Galactic winds in less massive galaxies are therefore not 
necessary to explain the MZ-relation. 

From the observational point of view, it is not clear if the 
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effective yield decreases or increases with galactic mass (see 
Tremonti et al. 2004; Erb et al. 2006a). New data on the redshift 
evolution of the gas fractions and of the effective yields have 
recently become available (Calura et al. 2008b; Mannucci et 
al. 2009). A forthcoming paper will be focused on the study of 
the redshift evolution of both the gas fractions and of effective 
yields, considering galaxies of various morphological types. 
2) At z = 0.07, the slope and the zero point of the observed 
MZ relation are affected by uncertainties related to the use of 
different metallicity calibration methods. Our predictions, ob- 
tained with models able to reproduce a large set of independent 
chemical evolution constraints for spirals and irregulars, may 
be useful to constrain the best calibration method. By means 
of our predictions, we suggest the best calibrations methods 
to be used by observers to derive gas metalhcities. Our results 
are consistent with the MZ relations obtained by means of the 
calibration of Kewley & Dopita (2002), Denicolo et al. (2002) 
and Pettini & Pagel (2004). 

3) From the redshift evolution of the MZ relation, we have 
seen how galaxies of various morphological types occupy 
different regions in the MZ plots at various redshifts. However, 
in these plots, in some cases, predictions for different galaxies 
tend to overlap. For this reason, the MZ plot alone cannot be 
used to constrain the morphology of the observed galaxies. 

4) To constrain the morphological type of the observed 
galaxies, a very helpful method is to consider, beside the MZ 
plot, O/H vs SFR and SFR vs mass diagrams. We have tested 
how the adoption of different values for the parameters used in 
this study, Zf and A,, influences the predictions at various red- 
shifts. At z = 3.5, z - 2.2 and z - 0.7, our models can repro- 
duce both the metallicities and the SFRs of the observed galaxies 
by assuming for all galaxies an average redshift of formation of 
Zf - 3. An age dispersion of A, = 3 Gyr for ellipticals, com- 
patible with the values coming from photometric studies of the 
Colour-magnitude relations of early-type galaxies in clusters and 
in the field. A; = 5 Gyr for spirals and A; = 10 irregulars, allows 
us to reproduce most of the observational constraints considered 
in this work. 

5) The contribution of each type of star forming galaxy to 
the MZ relation is a fuction of redshift. At z = 3.5, the observed 
stellar masses, metallicities and SFRs are reproduced by our 
models for eUiptical galaxies. At z = 2.2, the observed metallic- 
ities and SFRs indicate that the galaxies are likely to represent 
a morphological mix, composed partly by spirals (or proto- 
spirals) and partly by (proto-)eUipticals. At z = 0.7, we sUghtly 
underestimate the observed metalhcities. Furthermore, we tend 
to overestimate the observed SFRs, which, plotted against the 
stellar mass M«, show a weaker correlation than our predictions. 
At z = 0.07, we reproduce the shape of the observed MZ 
relation. However, the dispersion of the data is larger than what 
our predictions indicate. The discrepancies between predictions 
and observations at z = 0.7 and z - 0.07 could be due to the 
fact that, in our chemical evolution models, the star formation 
histories do not take into account stochastic events such as 
episodic starbursts. Such events could increase the dispersion 
of the observed metallicities and make the real star formation 
histories more irregular than the predicted ones. In any case, this 
does not represent a major concern since, as shown in Kewley & 
Elhson (2008), also the dispersion depends on the cahbration. 
Our predicted SFRs are compatible with the observed values of 
the SDSS galaxies, although this sample includes also objects 
with very high SFRs, up to log(S FR/Mo/yr) = 1.8 - 2, which 



we can not reproduce and which are sometimes even higher 
than the values observed at z = 0.7 in the GEMINI Deep Deep 
Survey (Juneau et al. 2005). A possible reason for such high 
SFRs could be due to overestimated dust attenuation in the 
SDSS sample (J. Brinchmann, private communication). 

6) At z = 3.5, the majority of the predicted metalhcity 
values for the ellipticals are higher than the observed ones by 
~ 0.5 dex. This difference is larger than the offset due to the 
use of different calibrations. By computing the dust surface 
mass densities for our elliptical models, we show that, if the 
AMAZE sample is composed mainly by the progenitors of 
local ellipticals, most of these galaxies may remain observable 
for times < 0.2 Gyr, considering also the uncertainties in the 
observed SFRs. After 0.2 Gyr, the bulk of galaxies may become 
heavily dust-obscured and hardly detectable by current surveys 
in the optical-UV bands. 

7) Our study of the MZ relation for various elements in 
spirals shows that the predicted dispersion at the lowest stellar 
masses depends on the lifetimes of the stellar progenitors 
producing the elements. The largest dispersions in the MZ plot 
are predicted for Si and Fe, which are produced by Type la SNe 
in considerable amounts. 

8) An analysis of the predicted [O/Fe] plots computed at 

z = 3.5 and at z = 2.2 for ellipticals, spirals and irregulars shows 
that, if interstellar Fe abundances in high-redshift star-forming 
galaxies were measurable, thanks to our chemical evolution 
models, it would be possible to have important hints on the 
morphology of these objects. 

Unfortunately, interstellar Fe abundance measures from 
emission lines are not achievable by means of present-day 
instruments. However, the determination of the iron abundance 
of star-forming galaxies at high redshift, as performed by 
Halliday et al. (2008), which use flie Rix et al. (2008) Fe 
absorption line, could represent an important result in this 
regard. The combination of these abundances with measures 
of the interstellar O abundance would provide the observables 
required by our study. 

It is also worth to note that Fe abundances can be obtained 
in absorption in gas-rich galaxies observed in the spectra of 
high-redshift QSOs, such as Damped-Lyman Alpha galaxies, 
together with the abundances of various other elements. A 
seminal work of Ledoux et al. (2006) points towards the 
existence of a possible mass-metallicity relation also in DLAs. 
In the future, this issue, which needs to be confirmed by further 
observations, wiU be an interesting subject for forthcoming 
theoretical chemical evolution studies. 
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Fig. 6. Predicted MZ-relation at z=0 for spirals and irregulars. In the figure is shown also the best fit to the local MZ-relation derived 
by Maiolino et al. (2008) from the data of Kewley & Ellison (2008). The open squares are the data by Lee et al. (2006), who 
measured the MZ relation in a sample of local dwarf galaxies. 
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Fig. 7. Predicted Mass-Metallicity for spiral galaxies (red contours) and iiTegular galaxies (blue contours) at z = 0.07. The solid 
lines of different colours are best-fit relations calculated by Kewley & Ellison (2008) using different metallicity calibrations. KK04: 
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Fig. 8. Redshift evolution of the MZ relation as a function of morphological type. The black points are the observations at various 
redshifts. The blue and red contours repesent the regions where 100 % and 75 % of our predictions lie, respectively. In the first, 
second and third column the results for ellipticals, spirals and irregulars are shown, respectively. From top to bottom and for each 
morphological type, predictions and observations at four different redshifts are shown: z - 0.07, z = 0.7, z = 2.2, z - 3.5. In this 
case, we have assumed a redshift of formation z/ = 3 and an age dispersion of A, = 3 Gyr for all galaxies. Small crosses at z = 0.07 
are from Kewley & Ellison (2008); diamonds and open triangles at z = 0.7 from Savaglio et al. (2005); the open squares at z = 2.2 
are from Erb et al. (2006a); at z = 3.5, the solid circles are from Maiolino et al. (2008). All observational data have been recalibrated 
as discussed in Maiolino et al. (2008). 
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Fig. 9. Redshift evolution of the MZ relation as a function of morphological type. Symbols and contours as in Fig. [8] In this case, 
we have assumed a redshift of formation zy = 6 and an age dispersion of A, = 3 Gyr. 
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Fig. 10. Predicted redshift evolution of the MZ, (O/H) vs SFR and SFR vs mass plots for ellipticals (green contours), spirals (red 
contours) and irregulars (blue contours) and as observed by various authors at z = 0.07, z - 0.7, z - 2.2, z - 3.5. For the predictions, 
we have assumed a redshift of formation Z/ = 3 and an age dispersion of A, = 3 Gyr for ellipticals. A, = 5 Gyr for spirals and 
A, = 10 Gyr for irregulars. Crosses at z = 0.07 from Kewley & Ellison (2008) for the M-Z plot, Kewley & Ellison (2008) and 
Brinchmann et al. (2004) for the (O/H) vs SFR plot, and Brinchmann et al. (2004) for the SFR vs mass plot. Diamonds and open 
triangles at z = 0.7 from Savaglio et al. (2005) and Maier et al. (2005). At this redshift, the SFRs are from Juneau et al. (2005). At 
z = 2.2, the open squares are from Erb et al. (2006a) and Erb et al. (2006b). At z = 3.5, the solid circles are from Maiolino et al. 
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Fig. 11. Upper panel: SFR vs stellar mass plot at z = 3.5. Green, magenta and cyan regions as above. The three dashed lines are 
isochrones at 0.05 Gyr, 0.1 Gyr, and 0.2 Gyr, from the left to the right of the plot. The three solid lines are the evolutionary tracks 
for the three baseline models for ellipticals considered in this paper. 

Lower panel: MZ plot at z = 3.5 as predicted for elliptical galaxies and as observed by Maiolino et al. (2008). The green, cyan 
and magenta regions represents our predictions for ellipticals with dust surface mass densities in the ranges cr^/ > Q.3Mq/ pc^, 
0.1 < cr^/0.3Mo/ pc^ < 0.3, cTrf < Q.IMq/ pc^, respectively. The three dashed lines are isochrones at 0.05 Gyr, 0.1 Gyr, and 0.2 Gyr, 
from the lowest line to the highest one. The three solid lines are the evolutionary tracks for the three baseline models for ellipticals 
considered in this paper. 
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Fig. 12. Predicted MZ relations for various elements in spiral galaxies. Here we have assumed zj - 3 and A, = 10 Gyr. 
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Fig. 13. Predicted [O/Fe] vs [Fe/H] at z = 3.5 (lower panel) and at z = 2.2 (upper panel) for elliptical galaxies (blue areas), spirals 
(cyan areas) and iiTegulars (magenta areas). Here we have assumed z/ = 3 and A, = 10 Gyr for spirals and irregulars and z/ = 3 and 
A, = 3 for ellipticals. 



